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Abstract 

The low-density expansion for a homogeneous interacting Bose gas at zero temper- 
ature can be formulated as an expansion in powers of \J pa?, where p is the number 
density and a is the S-wave scattering length. Logarithms of pa 3 appear in the co- 
efficients of the expansion. We show that these logarithms are determined by the 
renormalization properties of the effective field theory that describes the scattering of 
atoms at zero density. The leading logarithm is determined by the renormalization of 
the pointlike 3^3 scattering amplitude. 



The successful achievement of Bose-Einstein condensation of atomic gases in magnetic 
traps has created an explosion of interest in Bose gases of atoms. While a qualitative 
description of the condensation can be obtained using mean field methods 0, a more quanti- 
tative treatment requires including corrections from fluctuations around the mean field. The 
relative magnitude of these corrections grows with the number density of the atoms. They 
will therefore become more important as higher trap densities are achieved. 

In order to develop a deeper understanding of the fluctuations, it is worthwhile to go back 
to the simpler problem of a homogeneous gas of interacting bosons at zero temperature. This 
problem was studied intensively in the 1950's ||, ffl. A simple review was given by Yang in 
1960 0. The properties of the system can be calculated as an expansion in powers of \/ pa 6 , 
where p is the number density of atoms and a is their S-wave scattering length. For example, 
the expansion for the energy density has the form 

|(4tt-3V3) \og(pa 3 ) + K 
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where we have set Ti = 1. The coefficient of \] pd A was first obtained by Lee and Yang for 
a hard sphere gas ||. The coefficient of pa 3 \og(pa 3 ) was calculated by Wu, by Hugenholtz 
and Pines, and by Sawada Q. The correction npa 3 is the first term in the expansion that is 
sensitive to atomic parameters other than the scattering length. We have recently succeeded 
in calculating the constant k 0. 

In this Letter, we use a minimal subtraction renormalization scheme to deduce the general 
structure of the low-density expansion. We show that logarithms of pa 3 are related to the 
renormalization of the effective field theory that describes the scattering of atoms in the 
vacuum. In particular, the log(pa 3 ) term in (P is related to the renormalization of the 
3^3 scattering amplitude. We reproduce the leading logarithms in previous calculations 
using simple renormalization group methods. Our approach can also be used to determine 
the logarithms that appear at higher orders in the low-density expansion. 

Our starting point is an effective field theory that describes atoms with momenta 
much lower than their inverse size, which is on the order of the Bohr radius qq. Since the 



range of the interaction potential between 2 or more atoms is also on the order of a , the 
interactions appear pointlike on the scale of the de Broglie wavelengths of the atoms. The 
atoms can therefore be described by a field theory with a hamiltonian density that is a local 
function of the field: 

H = -^ f V 2 ^ + ~g&W + iflityty) 8 + (2) 
2m 4 3o 

For simplicity, we have assumed that the atoms have spin 0, so that they can be represented 

by a single complex field ip. The (ip^ip) 2 term represents 2 — > 2 scattering through an S-wave 

interaction with scattering length a given by g = 87ra/m, while the term (ifj^ip) 3 represents 

3^3 scattering. By adding additional terms that are higher order in the derivatives or 

in the number of fields, one can describe n — > n scattering of atoms in the vacuum with 

whatever accuracy is desired. In principle, the coefficients of these terms can be calculated 

from the n-body potentials that describe interatomic interactions. In the absence of such 

calculations, they can be taken as phenomenological parameters. 

By treating the interaction terms as perturbations, we can calculate the amplitudes for 

scattering of atoms with momenta on the order of p as an expansion in powers of pa$. 

This expansion is complicated by the presence of ultraviolet divergences. For example, the 

amplitude for the scattering of two atoms with momenta pi and P2, including the first 

perturbative correction, is 

mg r d 3 k 1 
9 [ ~ ~Y J (2tt) 3 k 2 - ( Pl + p 2 ) ■ k + Pl ■ p 2 - ie 

The integral, which is ultraviolet divergent, can be regularized by imposing a cutoff |k| < 

A. The linear divergence can then be cancelled by adding a counterterm proportional to 

mg 2 A(ip^if)) 2 to the effective hamiltonian (Q). The resulting expression for the scattering 

amplitude is rather complicated, as it includes terms that are suppressed by powers of pi/A 

and P2/A. A simple analytic result is obtained only in the limit A — » 00. 

A power ultraviolet divergence, such as the linear divergence in (||), indicates extreme 

sensitivity to short- distance atomic physics that is not accurately described by the effective 
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hamiltonian (|2]). A simple momentum cutoff is not an accurate model for the way atomic 
physics cuts off the momentum integrals. There is an alternative cutoff procedure, called 
"minimal subtraction", which is no more accurate a model for the cutoff, but provides an 
equally accurate description of the long-distance physics and has the virtue of simplicity. In 
minimal subtraction, linear, quadratic, and other power ultraviolet divergences are removed 
as part of the regularization scheme by subtracting the appropriate power of k from the 
momentum space integrand. In the case of the amplitude (H), 1/k 2 is subtracted from the 
integrand. The justification for this procedure is that the terms that are subtracted are 
dominated by short distances and can be cancelled by counterterms in the Hamiltonian. In 
minimal subtraction, logarithmic ultraviolet divergences are treated differently from power 
divergences. This is reasonable, because logarithmic ultraviolet divergences represent real 
physical effects, while power ultraviolet divergences are simply artifacts of the regularization 
procedure. This difference is reflected in the fact that the coefficient of a power divergence 
A p depends on the regularization prescription, while the coefficient of log(A) does not. The 
reason for this is that the logarithm of A must match onto the logarithm of some physical 
momentum scale, and therefore its coefficient has a real physical meaning. We regularize 
logarithmic ultraviolet divergences by imposing a cutoff |k| < A on loop integrals. After using 
renormalization to remove divergences from subdiagrams, we isolate the divergent terms 
proportional to log(A) and then take the limit A — > oo in the remainder. The cutoff A is called 
the "renormalization scale" . With minimal subtraction, all power divergences and those parts 
of logarithmic divergences that arise from momenta greater than the renormalization scale 
A are absorbed into the coupling constants in the effective hamiltonian. 

The advantage of minimal subtraction is that it makes it much easier to disentangle the 
effects of different momentum scales in multiloop diagrams. With a conventional momen- 
tum cutoff, a diagram can be an extremely complicated function of the cutoff A, the external 
momenta, and the momentum scales that can be formed from the parameters in the hamil- 
tonian. With minimal subtraction, the possible dependence of a diagram on the cutoff is 
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greatly simplified. The dependence can only be polynomial in log(A), with the logarithms 
arising from logarithmically ultraviolet divergent subdiagrams. This makes it much easier to 
analyze the divergences in a multiloop diagram. The relative simplicity of minimal subtrac- 
tion is illustrated by the fact that it gives a simple expression for the scattering amplitude (^) 
that is independent of the renormalization scale A: g[l+img\pi — P2 1 /(167r)] . The renormal- 
ized parameter g is independent of A, and satisfies a trivial renormalization group equation: 
A(d/dA)g = 0. With a conventional momentum cutoff A, the amplitude (^) is a complicated 
function of pi, P2, and A. For A 3> |pi|, |p 2 |, it reduces to g(A)[l — img\pi — P2 1 / (87r)] , 
where g{A) is a renormalized coupling constant that satisfies A(d/dA)g(A) = —mg 2 A/(47r 2 ). 
The running of g{A) is generated by a power ultraviolet divergence and therefore has no 
real physical significance. The scale-invariant parameter g defined by minimal subtraction 
provides an equally accurate description of the long-distance physics. 

In the vacuum, the simplest quantity in which logarithmic ultraviolet divergences ap- 
pear is the 3^3 scattering amplitude. There is a tree-level contribution from the (if)*ijj) 3 
term in the effective hamiltonian, but there are also additional contributions that involve 
successive 2^2 scatterings. They include the 2-loop diagrams shown in Figure 1, which 
involve 4 successive 2^2 scatterings. These diagrams are logarithmically ultraviolet diver- 
gent. Removing the linear ultraviolet divergence from a sub diagram of the first diagram in 
Figure 1 by a subtraction in the integrand and then imposing a cutoff A, we find that the 
logarithmically divergent term is — 3(47r — 3\^3)m 3 g 4 log(A)/(327r 3 ). 

The renormalization scale A represents an arbitrary separation between short-distance 
effects, which are taken into account through the parameters in the effective hamiltonian 
and long-distance effects, which are calculated using the effective theory. Physical quantities, 
such as the 3^3 scattering amplitude, should therefore be completely independent of 
A. The explicit A-dependence from the two- loop diagrams must therefore be cancelled by 
implicit A-dependence from the coefficient #3 in the tree-level contribution. This statement 
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can be expressed as a "renormalization group equation" : 

A ^ 3(A) = ^ 47r - 3v/ 3) m ^ 4 • ^ 

It tells us that the parameter g 3 is really a "running coupling constant" that increases 
logarithmically as the momentum scale is increased. The renormalization scale A can be 
interpreted as the inverse of the spacial resolution. As this resolution is decreased, we resolve 
part of the "pointlike" 3 — > 3 scattering amplitude into the successive 2 — > 2 scatterings 
represented by the diagrams in Figure 1. The contributions from the two diagrams have 
opposite signs and the net effect is that the coupling constant g% decreases as A decreases. 

If the running coupling constant g% is determined at the scale 1/ao of atomic structure, 
it can be calculated at a lower momentum scale A by solving the renormalization group 
equation ([J): 

2 3 (A) = g 3 (l/a ) - -|^(47r-3v / 3)mVlog(^) . (5) 

Regardless of the sign of g^(l/ao), (73(A) eventually turns negative for sufficiently small A. In 
describing the scattering of atoms with momenta on the order of p, it is appropriate to choose 
the renormalization scale A to be of order p. The coefficients in the perturbation expansion 
can include logarithms of the form log(A/p), which are generated by logarithmically divergent 
subdiagrams. By choosing A to be of order p, such large logarithms are removed from the 
coefficients and absorbed into the parameters of the effective hamiltonian. 

We now consider the energy density £ of the Bose gas with number density p. In order 
for the system to have a homogeneous ground state that is stable, or at least metastable, the 
scattering length a must be positive. Neglecting for the moment the effects of fluctuations, 
the field ip develops a vacuum expectation value yfp. The energy density at tree level is 

So = \ 9p2 + 3^ 3(A)p3 + ' • ' • (6) 

Setting g = 8iTa/m, where a is the S-wave scattering length, the first term above reproduces 
the leading term in ([[]). The 1-loop contribution is the sum of the zero-point energies of 
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the Bogoliubov modes, with power ultraviolet divergences removed by subtractions in the 
integrand: 



where e(k) = kyfk 2 + 2mgp/ (2m). This integral reproduces the first correction term in ([!]). 

The correction of order pa 3 in (fj) requires the calculation of two-loop diagrams. However 
the term proportional to pa 3 log(pa 3 ) can be obtained without any further calculation. The 
reason is that this term is related to the renormalization of the amplitude for 3 — > 3 scattering 
in the vacuum. The two-loop diagrams for the energy density contain logarithmic ultraviolet 
divergences. From the expression for the Bogoliubov energy, we see that the momentum 
scale associated with the quasiparticles modes is \J1mgp. The logarithmic ultraviolet diver- 
gences from the two-loop diagrams will therefore be proportional to \og(A/ ^2mgp). Large 
logarithms such as this in the coefficients in the perturbation expansion can be avoided by 
choosing the renormalization scale A to be on the order of y / 2mgp. With this choice of the 
renormalization scale, all such logarithms are absorbed into the parameters in the effective 
hamiltonian. Substituting A = y/Wnap in (||) and using the expression for g 3 (A) in (|5|), 
we reproduce the term containing the logarithm in (JlJ) . We have determined the constant 
k under the logarithm by calculating the 2-loop Feynman diagrams for the energy density 
explicitly. The details of the calculation will be reported elsewhere ||. As noted previously 
Q|, the pa 3 term is the first term in the low density expansion ([J) that is sensitive to atomic 
physics parameters other than the scattering length a. The only additional parameter that 
enters at this order is the pointlike 3 — > 3 scattering amplitude g%. 

The renormalization group together with minimal subtraction can also be used to deter- 
mine the leading logarithms in the low density expansions for other quantities. Corrections 
to the sound velocity have been calculated by Beliaev || , including the V pa 3 term and the 
pa 3 log(pa 3 ) term. The logarithm can be obtained by calculating the logarithmic ultraviolet 
divergences in the two-loop corrections to the propagator. Alternatively, it can be obtained 
trivially using the methods described above. Including the correction from the (ip^ip) 3 term 
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in (0), the sound velocity at tree level is v 2 = gp/(2m) + g 3 (A)p 2 /3. Setting A = ^16nap 
and using the expression (|5]) for (73(A), we reproduce the pa 3 log(pa 3 ) correction calculated 
by Beliaev. 

These methods can also be used to determine the coefficients of the logarithms that appear 
at higher orders in the low-density expansion for the energy density. For example, there is a 
(pa 3 ) 3 / 2 log(pa 3 ) correction to the energy density which arises from logarithmically divergent 
two-loop subdiagrams in 3-loop diagrams. This terms can be determined easily by taking 
into account the (ip^ip) 3 term in the Bogoliubov energy: e 2 (k) = k 2 {k 2 + 2mg p) / (4m 2 ) + 
g^(A)p 2 k 2 /3. Inserting this into the 1-loop expression (|7]) for the energy density, expanding to 
first order in g 3 , and using the expression (|j) for g 3 (^167iap), we obtain the (pa 3 ) 3 / 2 log(pa 3 ) 
correction. 

Thus far we have only considered logarithms in the low density expansion that are related 
to the renormalization of the 3 — > 3 scattering amplitude. At higher orders in the low density 
expansion, there are also logarithms that are related to the renormalization of other terms in 
the effective lagrangian, such as the term g^t/j^if)) 4 which describes 4-^4 scattering through 
a point-like interaction. The 4 — » 4 scattering amplitude includes logarithmically ultraviolet 
divergent corrections from 4-loop diagrams that involve 7 successive 2^2 scatterings and 
also from 2-loop diagrams that involve three 2^2 scatterings and a 3 —>■ 3 scattering. 
The explicit A-dependence from these loop diagrams must be cancelled by the implicit 
A-dependence from the pointlike 4-^4 scattering amplitude (74(A). As a consequence, 
the renormalization group equation analogous to (^) for (Arf/rfA)^ includes terms on the 
right side that are proportional to m 6 g 7 and m 3 g 3 g 3 . The solution for 54(A) analogous to 
(||) includes a term proportional to m 6 g 7 log 2 (Aao). Choosing A = ^/IQ-nap to avoid large 
logarithms from loop diagrams, we find that the term ^p 4 in the mean-field expression for 
the energy density gives rise to a correction to ([!]) that is proportional to (pa 3 ) 2 log 2 (pa 3 ). 

A renormalization group analysis of the dilute Bose gas at nonzero temperature was 
recently carried out in Ref . || . The authors derived renormalization group equations for the 
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chemical potential p and for the coupling constant Vq = 2g using a conventional momentum 
cutoff. If they had used minimal subtraction, their analysis would not have been modified 
dramatically. Their equations for dp/dl would remain unchanged. In their equations for 
dVo/dl, there would be an additional term V^ 2 mA/(27r 2 ) on the right side which cancels the 
leading power of A in the equation at T = 0. The renormalization group trajectories for these 
two renormalization schemes would differ significantly only near the initial cutoff, where both 
calculations would be dominated by cutoff artifacts. They would be essentially identical near 
the critical point for Bose condensation. The logarithmic evolution of the coupling constant 
(73 was not seen in the analysis of Ref. ||, because they considered renormalization effects 
from one-loop diagrams only. 

In this Letter, we have shown how the structure of the low-density expansion for a 
Bose gas is determined by the renormalization properties of the effective hamiltonian that 
describes the scattering of atoms at zero density. The low density expansion for the energy 
density has the general form 

e = -EE c «i (pa 3 r /2 log V 3 ) • (8) 

m n=0l=0 

The maximum power of the logarithm has been determined to be l n = 0, 0, 1, 1, 2 for n = 
0, 1, 2, 3, 4, respectively. The dimensionless coefficients C n \ are polynomials in the generalized 
coupling constants of higher order terms in the effective hamiltonian (Q), with only a finite 
number of these coupling constants appearing at any given order in \j pa 6 . The coupling 
constant first appears at order pa 3 . Additional coupling constants enter at order (pa 3 ) 3 / 2 . 
The general structure in (|8|) follows automatically from the renormalization group together 
with the minimal subtraction renormalization scheme. This powerful method should also be 
useful for analyzing the corrections from fluctuations around the mean field for atomic gases 
in magnetic traps. 

This work was supported in part by the U. S. Department of Energy, Division of High 
Energy Physics, under Grant DE-FG02-91-ER40690. 
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Figure Caption 

Figure 1. Two- loop diagrams that give logarithmically ultraviolet divergent contributions to 
the 3^3 scattering of atoms. 
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